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;h ' Abstract 



We discuss the properties of hadronic fluctuations, i.e. fluctuations of net quark 
and isospin numbers as well as electric charge, in the vicinity of the QCD transition 
in isospin-symmetric matter at vanishing quark chemical potential. We analyse 
second- and fourth-order cumulants of these fluctuations and argue that the ratio of 
quartic and quadratic cumulants reflects the relevant degrees of freedom that carry 
the quark number, isospin or charge, respectively. In the hadronic phase we find 
that an enhancement of charge fluctuations arises from contributions of doubly 
charged hadrons to the thermodynamics. The rapid suppression of fluctuations 
seen in the high-temperature phase suggests that in the QGP, net quark number 
and electric charge are predominantly carried by quasi-particles, with the quantum 
numbers of quarks. 
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1 Introduction 



In a medium of strongly interacting elementary particles, described by equi- 
librium thermodynamics of Quantum Chromodynamics (QCD), the tran- 
sition between the low-temperature hadronic and high-temperature quark- 
gluon phase occurs in a small temperature interval and leads to a rapid change 
in entropy and energy density. Although it is, most likely, that this transi- 
tion is not a true phase transition, but rather a smooth but rapid crossover, 
it is expected that this transition leads to large fluctuations of energy and 
net quark number densities. In fact, lattice studies of various hadronic sus- 
ceptibilities, such as net quark number, isospin or charge susceptibilities, 
show that quadratic fluctuations of these quantum numbers increase rapidly 
in the transition region [1, 2]. These susceptibilities continue to grow also 
above the transition temperature, Tc, and come close to the ideal-gas value at 
temperatures T ~ 2Tc. Below Tc bulk thermodynamics as well as the generic 
structure of susceptibilities is found to be in good agreement with properties 
of a hadronic resonance gas [3]. 

Recent lattice studies of 2-flavour QCD at non-zero quark (/Xg) and isospin 
(/i/) chemical potential [4-7] have sustain that higher-order derivatives of the 
QCD partition function (generalized susceptibilities), show even more pro- 
nounced variations with temperature in the vicinity of the QCD transition 
temperature. In particular, the fourth-order derivatives with respect to Hqj, 
or equivalently with respect to the corresponding up and down quark chemi- 
cal potentials (/i^, l^d), have pronounced peaks at the transition temperature. 
As a consequence quark number fluctuations calculated at non-zero /ig in- 
crease with increasing values of fig and develop a cusp along the transition 
line. This is in agreement with expectations based on the (mean-field) anal- 
ysis of quark number fluctuations in sigma models [8] as well as NJL models 
[9]. These models share the relevant 0(4) symmetry with 2-flavour QCD, and 
give some hint for the existence of a second-order phase transition point in 
the QCD phase diagram at some non-zero vahic of the quark chemical poten- 
tial. In this paper we will, however, concentrate on a discussion of hadronic 
fluctuations at vanishing quark chemical potential. This is a regime of the 
QCD phase diagram most relevant to current heavy-ion experiments at RHIC 
as well as future experiments at the LHC. 

Lattice calculations of bulk thermodynamics, e.g. of the equation of state, 
as well as studies of generalized susceptibilities, show signiflcant deviations 
from the ideal-gas behaviour for a rather large temperature range above 
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Tc- In particular, the equation of state shows sizeable deviations from the 
ideal-gas relation, e = 3p, up to T ~ (3-4)Tc. This is in line with recent 
experimental findings that hint at the formation of a strongly interacting 
medium in heavy- ion collisions at RHIC energies [10]. 

It has been recently argued that this strongly interacting medium may, 
in fact, consist of a large set of heavy coloured bound states, which could 
exist in a temperature interval Tf. < T<(l-2)Tc [11, 12]. The so-called sQGP 
model suggests that these bound states can contribute as much as 20% to 
the pressure in the high-temperature phase of QCD. Such states clearly will 
also contribute to hadronic fluctuations. We will discuss here to what extent 
fluctuations of such states are consistent with lattice results on quark number 
and charge fluctuations in the QGP. 

In the next section we will give a definition of the observables we are 
going to analyse and present the lattice results we are going to discuss in 
the following sections. In section 3 we discuss the properties of quadratic 
and quartic fluctuations in the vicinity of the QCD transition temperature. 
Section 4 is devoted to a discussion of hadronic fluctuations in the low- and 
high-temperature phases of QCD. Section 5 contains our conclusions. 

2 Hadronic fluctuations at /x^ = 

Quark number as well as charge and isospin fluctuations are obtained from 
the grand canonical partition function of QCD through appropriate combi- 
nations of derivatives with respect to quark chemical potentials. As we will 
concentrate here on a discussion of their properties at finite temperature and 
vanishing quark chemical potential in 2- flavour QCD, we will use the QCD 
partition function with non-vanishing up and down quark chemical poten- 
tials only as a tool to generate the derivatives of interest. After having done 
so the chemical potentials will be set to zero. 

The pressure of a thermodynamic system is directly related to the grand 
canonical partition function: 

p(y,r,/x„,//d) = ^inZ(y,r,//„,/Xd) , (i) 

where /x„ and iid are the chemical potentials for up and down quarks, respec- 
tively. 

The fluctuations of quark number, isospin and charge as well as expec- 
tation values of higher-moments of these observables can be obtained from 
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Figure 1: Quadratic fluctuations of tlie quark (left) and isospin numbers. 



generalized susceptibilities: [6, 7, 13] 

an+m p 

n,m ^ _^ f_ _ /2) 

These susceptibilities, evaluated for vanishing quark chemical potentials, de- 
fine the coefficients of a Taylor series for the pressure expanded simultane- 
ously in terms of [Xu and [Xd- In the following we will consider only the mass- 
degenerate case (m„ = nid) for which the susceptibilities are symmetric, 
^n,m _ ^m,n_ f^g^gnt lattice calculations have shown that the fourth deriva- 
tives, as well as the off-diagonal susceptibilities, x^'^, become large in the 
vicinity of the transition temperature [6]. Quartic fluctuations thus seem to 
be particularly sensitive to details of the transition to the high-temperature 
phase. 

Fluctuations of net quark {Nq = AN^ + ANd), isospin (Nj = ANu-ANa) 
and charge {Q = ^ANu — ^ANd) numbers, with AA^^^^ = N^^—Nu.d denoting 
the difference of quarks and antiquarks, are conveniently discussed in terms 
of quark and isospin chemical potentials: 

/^g = — ^ — , = — ^ — • (3) 

Charge fluctuations can then be obtained from appropriate combinations of 
derivatives with respect to ^u,d or /^g,/, respectively: 

d 2d Id Id 
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The first two, non- vanishing derivatives with respect to H:q,i,Q then yield 
,^ 1 9^1nZ 1 ,2\ T ^ 

■'J^FfSag^ = ^«(^iv.r>-3((wjV) . (6) 

We note that the derivatives of the partition function with respect to //g 
are directly related to the expansion coefficients calculated in Refs. [5, 6]. 
Furthermore, related to the expansion coefficient Cg, while can be 
obtained from a combination of C2 and C2, which were obtained also in these 
references. 

As is evident from Eq. (4) charge fluctuations can easily be expressed 
in terms of corresponding fluctuations in quark and isospin (or net up and 
down quark) numbers. While the quadratic charge fluctuations arc entirely 
determined by the corresponding quark number and isospin fluctuations, the 
quartic charge fluctuations are also sensitive to correlations between 5Nq and 
5Nj. They thus provide independent information: 

(f, = l(dl + 8l4 + 54(((57V^)2(57V,)2)-3W7V,)2)) . (6) 

In Fig. 1 we show the two independent second-order cumulants d\ and (ig- 
They are directly proportional to the coefficients C2 and C2 calculated in 
Refs. [5, 6]. The three fourth-order cumulants^ are shown in Fig. 2. 

While the second-order cumulants show the well known continuous in- 
crease towards the high-temperature ideal-gas value, the fourth-order cu- 
mulants clearly have a pronounced peak at the transition temperature and 
exceed the ideal-gas value by more than a factor 4. Here it is worthwhile to 
note that the lattice results shown in Figs. 1 and 2 are based on calculations 
that were performed with still rather heavy quarks, corresponding to a pion 
mass of about 770 McV. Although we still no have satisfactory understand- 
ing of the universal behaviour of the QCD transition in 2-flavour QCD, the 
studies of the quark mass dependence of, for instance, the chiral susceptibil- 
ity [14-17] suggest that the peaks in susceptibihties increase with decreasing 

^ These results are based on fourth-order derivatives of the pressure calculated in 
Refs. [5, 6]. 
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Figure 2: Fourth cumulant of the quark (top, left) and isospin (top, right) 
number fluctuations as well as the fourth cumulant of the charge fluctuations 
(bottom, left). 



quark mass. We thus also expect that the peaks found in fourth derivatives 
of the pressure will increase with decreasing quark mass. In the next section 
we briefly discuss the generic structure of second- and fourth-order cumulants 
expected to be found in the chiral limit. 
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3 Chiral symmetry restoration and 
hadronic fluctuations 



In the vicinity of the chiral transition temperature, and for small values of the 
quark chemical potentials, the pressure receives contributions from a regular 
and a singular part: 

P{T, l^u: IJ-d) = Pr{T, Hu, /^d) + Ps{i, fiu: fid) , (7) 

where we have introduced in the singular part the reduced temperature 
i = \T — Tc\/Tc and the dimensionless chemical potentials fiu,d = fJ'u,d/T. 
Deviations from criticality are controlled by a generalized version of the re- 
duced temperature variable, which takes into account that the temperature- 
like direction at non-vanishing quark chemical potential in general will be a 
mixture of the couplings t and /iq/T. This generalized reduced temperature 
has to respect the invariance of the QCD partition function under exchange 
of particles and antiparticles {/iq — > —/iq)'- 

t = t + A{jil + /i^) + Bfiufid 

= t + {2A + B)jJil + {2A-B)fij . (8) 

The singular part of the pressure may then be written as 

Psil ftu: fid) ~ b~'fs{b'/^H) ~ , (9) 

where we have introduced an arbitrary scale parameter b and where yt = 2— a 
is the usual thermal critical exponent, which is related to the critical exponent 
a. Equation (9) shows that the critical behaviour of hadronic fluctuations 
= A*(i = is controlled by the critical exponent of the specific heat, a [8] . 
The contribution of the singular part to second- and fourth-order cumulants 
of quark {k = q) and isospin {k = /) is then given by 

(d'Ps/df^l) ~ t'-\ 

(d'p^idi^i) ^ , - t--- (10) 

For a conventional second-order phase transition with a critical exponent 
a > we would thus expect that the fourth-order cumulants diverge in 
the chiral limit. However, if the chiral transition is indeed of second-order 
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we expect it to belong to the universality class of three-dimensional 0(4) 
symmetric spin models. In this case the exponent a will be negative, a ~ 
—0.21 [18]. The contribution of the singular part will then vanish at t = and 
the value of the cumulants at t = will thus be determined by the regular 
part alone. In the vicinity of i = the temperature dependence of the 
second-order cumulant will then be dominated by the regular part. This is 
in accordance with the rapid but smooth increase of second cumulants across 
the transition temperature. The temperature dependence of the fourth-order 
cumulant, however, will be dominated by the singular part and as a is small 
the temperature variation will result in a sharp cusp. This is in accordance 
with the behaviour seen in Fig. 2. 

If the chiral transition in 2-flavour QCD is of second-order, we thus expect 
the second- and fourth-order cumulants to stay finite at the critical point. 
This, however, only holds for the critical point at \iq = 0, which we have 
discussed above. If the transition in 2-flavour QCD a,t /iq — indeed belongs 
to the 0(4) universality class, this transition point will be part of a line of 
second-order phase transitions at non-zero quark chemical potential. At a 
transition point with fig > 0, one would then expect that already the second- 
order cumulant of quark number fluctuations develops a cusp [8] and that 
the fourth-order cumulant would diverge. 

4 Hadronic fluctuations in a thermal medium 

In the previous section we have discussed the qualitative structure of second- 
and fourth-order cumulants in the vicinity of the QCD transition tempera- 
ture. Here we want to look in somewhat more detail into the temperature 
dependence of these cumulants in the low and high-temperature phases of 
QCD. In particular, we want to argue that the ratios of fourth- and second- 
order cumulants of quark number and charge fluctuations, 

are sensitive observables that allow the identification of the carriers of quark 
number and electric charge in a thermal medium. We will first discuss this 
for a hadron gas at low-temperature and then move on to a discussion of the 
relevant degrees of freedom in the high-temperature phase. This discussion 
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can easily be extended to the case of isospin fluctuations, which, however, 
we will not consider here any further. 

4.1 Fluctuations in hadronic matter 

We describe the low-temperature phase of QCD in terms of a resonance 
gas with free particle dispersion relations for all constituents [3, 19]. The 
resulting pressure, however, also accounts for interactions between hadrons, 
to the extent that the thermodynamics of an interacting system of elementary 
hadrons can effectively be approximated by that mixture of ideal-gases of 
stable particles and resonances [19, 20]. The pressure P is then expressed as 
a sum over all mesonic and baryonic degrees of freedom: 



Using Eq. (3) to replace the up and down quark chemical potentials by 
and /I I, the contribution of mesons, baryons and their resonances belonging 
to fixed baryon and isovector quantum numbers is obtained from Eq. (12) 
in compact form. With the exception of the pion, all hadrons are heavy, 
with respect to the temperatures of interest. For these particles quantum 
statistics plays no role, and the Boltzmann approximation of the partition 
functions is appropriate. This yields, in the meson sector. 



^;^=EF(r,m.)oosh?^ , (13) 

i 

where the sum is taken over all non-strange mesons^, being the third 
component of the isospin of the species i. The function F{T, rrii) is the 
contribution of particles and the corresponding antiparticles of mass rrii with 
spin degeneracy factor di to the pressure at vanishing baryon and isospin 
chemical potential 



F(r.„.,^|i^^(f)V,(^) . (14) 



^In order to be able to compare with lattice calculations performed for 2-flavour QCD, 
we have neglected the contribution of strange particles to the hadronic pressure. 
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Neglecting the mass difference between isospin partners, the meson con- 
tribution arises from isospin sing let (/f ^ = 0) and triplet G^^) (/f ) = 
0, ± 1) mesons, respectively. Thus, 



^""^^yf '^'^ - G^'Kt) + lG(^)(r)[2cosh(^) + 1] . (15) 

The contribution to the pressure arising from non-strange baryons and 
their resonances is given by 



^!^^^ = i:F(r,™.)cosh5^i±|^ . (16) 

i 

The baryonic part of the pressure receives contributions from isospin dou- 
blet F(2) ( /f^ = ±1/2) and quartet F^^) (/f^ = (±1/2, ±3/2)) baryons, 
respectively. Again ignoring the mass differences between isospin partners, 
the total baryonic contribution to the pressure is given by 



+ iF("(r)cosh(^)[cosh(^)+cosh(^)] . (17) 

We note that only the isospin quartet baryons contain doubly charged hadrons, 
which is reflected in the appearance of a term proportional to cosh(3////T). 

Within this resonance gas formulation we now can evaluate cumulants of 
the quark number, isospin and charge in the confined phase of QCD. In fact, 
by examining the cumulants introduced in section 2, we can find the con- 
tributions of different hadronic quantum number channels and can establish 
relations between the various cumulants. Within the Boltzmann approxi- 
mation for the hadron resonance gas, which we have introduced above, we 
find 



dl{T) = 8l(F(2)(r) + F(^)(r)) , (18) 
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d'iiT) = ^G(2)(r) + V)(T) + ^F(^)(T) , 

d^{T) = ^G'(')(T) + ^F(2)(r) + ^F(^)(r) . (19) 

(j Zi 

These equations are valid for isospin-symmetric systems at vanishing baryon 
number density. The corresponding derivatives of the contributions to the 
total pressure (Eqs. (13-17)) were thus taken at /x/ = /Xg = 0. The resonance 
gas model leads to simple relations for the ratios of cumulants of quark 
number and charge fluctuations. In particular, we found 

which is shown to agree well with the lattice results below Tc [5, 6]. 

The ratio i?4 2 is completely insensitive to details of the hadron spectrum 
contributing to the resonance gas. In particular, it reflects the fact that the 
quark number is carried only by baryons, i.e. the relevant degrees of freedom 
either carry net quark number 0, in which case they do not contribute to 
i?4 2, or 3, in which case they contribute to with a relative factor of 3^. 

This is similar for the ratio of cumulants of charge fluctuations, although 
we here have to take into account that hadrons with charge 0, 1 and 2 can 
contribute to the thermodynamics. The ratio become sensitive to the 
relative abundance of doubly charged hadrons: 

^'''-d^ " my) ^^^^^^^ 

4^(3) + 3F(^) + 27FW 
- 4^(3) + 3F(2) + 9FW ■ ^ 

At low-temperature, this ratio is dominated by charge fluctuations in the pion 
sector, which contributes to G^^\ The ratio will thus approach unity at low- 
temperatures and monotonically increase with temperature. This indicates 
that in the low-temperature limit all charged degrees of freedom carry one 
unit of charge. The ratio R^2 does increase only from contributions arising 
from isospin quartet baryons, which can carry charge Q = 2. In fact, denoting 
the contribution of all hadrons carrying charge Q = ±n by F^~"{T), one can 
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verify that Eq. (21) can be rewritten as^ 



Results from lattice calculations for the ratios defined in Eqs. (20) and 
(21), are shown in Fig. 3. In the low-temperature phase these ratios clearly 
show the basic features expected from the analysis of a hadron resonance 
gas. Above Tg, however, they rapidly drop and strongly deviate from the 
resonance gas behaviour. In fact, already at T ~ 1.5Tc the ratios are close to 
the corresponding ideal-gas values. In high-temperature perturbation theory 
the leading 0{g^) corrections to quadratic and quartic cumulants cancel in 
the corresponding ratios: 



^l2(^)|pe..th. - ^ + 0(^')' (23) 



RUt)^,..... - + (24) 



These values are also shown as solid lines in Fig. 3. 



4.2 Fluctuations in the QGP 

The considerations presented in the previous section for hadrons below T^, 
carry over to any thermodynamic system containing heavy quasi-particle 
states that are build up from quarks and gluons. As long as their masses 
are significantly larger than Tc, their contribution to the thermodynamics 
can be treated within the Boltzmann approximation. States with mass 
containing a net number of A;„ up quarks and down quarks will contribute 
to the pressure of the system with 

= F<(r, mi) cosh {kuliu + k^iid) , (25) 

where F{T,mi) is given by Eq. (14), i.e. it again combines the contribution 
of a particle and its antiparticle. 

"^This becomes obvious when one reconstructs the origin of the algebraic factors appear- 
ing in Eq. (19). For instance, the factor 3/2 in front of G^^^ indicates that, out of 3 triplet 
states, 2 carry charge IQj = 1 and the 3rd state docs not contribute as it carries charge 
(5 = 0. Similarly the factors 3/2 and 9/2 appearing in front of the quartet contribution 
F^^^ arise from the fact that, out of the 4 states, 2 carry charge \Q\ = 1, 1 state carries 
charge Q = 2 and the 4th state does not contribute as it carries charge Q = 0. 
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Figure 3: The ratios of fourth and second cumulants of quark number (left) 
and charge (right) fluctuations. At low-temperatures the ratios i?! 2 and i?^2 
are compared to resonance gas model predictions. At high-temperature the 
asymptotic ideal-gas value is indicated by a horizontal line. The dashed line 
represents the estimate for the value of the ratios at T = 1.4Tc in the presence 
of coloured qq and qg states asssuming that the qq states contribute with half 
their statistical weight as discussed in the text. 



A model containing heavy coloured bound states above has been put 
forward by E. Shuryak and I. Zahed [11] to explain the experimental evidence 
for strong interactions in the hot and dense medium of quarks and gluons 
created in heavy-ion collisions at RHIC. In the strongly coupled QGP model 
(sQGP), it is assumed that, in addition to heavy quasi-particles (quarks and 
gluons with large thermal masses), also a large number of coloured bound 
states contribute to the thermodynamics in the high-temperature phase of 
QCD. The model suggests, that aside from states carrying net quark number 

1 and charge Q = 1/3 and 2/3 there will also be states with net quark number 

2 and charge Q = \ and 4/3 contributing to the thermodynamics in the QGP. 
From the discussion given in the previous section it should be obvious that 
this must have consequences for quark number and charge fluctuations in the 
QGP and the ratios ^ and i?^2 should be sensitive to these new states. 
The lattice results shown in Fig. 3, can thus put bounds on their relative 
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abundance. 

For our discussion of baryon number and charge fluctuations, the quark- 
quark and quark-gluon bound states are of particular interest. The sQGP 
model suggests that the latter melt at Tmeit = Tqg ~ 2Tc while the former are 
only weakly bound and melt already at Tmeit = Tqq ~ lATc [11, 12]. It has 
been argued in Ref. [11] that the gradual disappearance of these states could 
be taken into account by reducing their contribution to the thermodynamic 
potential with increasing temperature. For the corresponding weight factor, 
we will use here the parameterization introduced in Ref. [11]: 

1 + exp (6(T - Tmeltj/Tc) 

where C controls the temperature interval over which the contribution of a 
given bound state to the partition function gradually disappears. 

In order to estimate the contribution of coloured states to the fluctuations 
of baryon number and charge, we consider the contribution of these states to 
the density-, (chemical potential)-dependent part of the pressure, 

^ = ^(F,(T) + i?2(T,2)F,,(T))(cosh(/x„/T) + cosh(/.,/T)) (27) 
+^Rc{T, lA)Fgg{T){cosh{2fXu/T) + cosh{2nd/T) 

+ cosh((//„ + /Xd)/r)) , 

where we again used the Boltzmann approximation as all states are expected 
to have masses > lOTc [11]. We will approximate the masses of quark- 
gluon and quark-quark bound states as twice the quasi-particle mass of 
quarks, and will thus neglect mass shifts arising from a non-vanishing bind- 
ing energy. Moreover, we used C = 2 for the quark-gluon bound states [11] 
and will here only discuss the contribution of quark-quark bound states as a 
function of Rc- With this we can estimate the contribution of quark-quark 
states to the quadratic and quartic quark number and charge fluctuations 
introduced in Eq. (5): 

dl = F,(r) + (i?2(r,2)F,,(r) + 16i?c(r,1.4)F,,(r)) . (28) 



and 



d^ = ^(^^.(r) + (^^2(T,2)F,,(T) + Hi?c(T,1.4)F,,(T))) 
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^ ^(yF,(T)+(^i?2(T,2)F,,(r) + i|^i?c(T,1.4)F,,(T))) (29) 

Apparently both ratios Rl'2 are enhanced over the corresponding ratios that 
would only contain states with the quantum number of quarks"^ owing to the 
presence of states that carry two units of quark number. Assuming that the 
coloured bound state masses are about twice as large as the quasi-particle 
quark masses, we find at the melting temperature of the quark-quark states 
RI2 — 1.55 and i?^2 = 0-52. These values are shown as the dashed lines 
in Fig. 3. How quickly this enhancement disappears at larger temperature 
depends on the precise choice of the suppression factor Re- At Tmcit — l-4T'c, 
the contribution of quark-quark states should be suppressed below 10% of 
that of an ordinary bound state {Rc ^ 0.1) in order to be compatible with 
the lattice results. This also could mean that the melting temperature is 
substantially smaller (Tmcit ^ l-lT'c). 

We also note that the perturbative values given in Eq. (24) are obtained 
for a massless quarks-antiquarks (fermion) gas. These ratios are smaller by 
about 50% than the corresponding values in the Boltzmann approximation. 
This indicates the uncertainties inherent in the estimates given above and 
suggests that if coloured bound states exist in the QGP and contribute to 
the thermodynamics, this presumably has to be discussed in terms of a more 
complex interaction pattern, which also takes into account quantum effects. 

5 Conclusions 

We have discussed some generic features of quark number, isospin and charge 
fluctuations in 2-flavour QCD. We argue that the ratio of quartic and quadratic 
fluctuations are sensitive observables that can directly provide information 
on the constituents of the thermal medium that carry net quark number and 
electric charge, respectively. We have shown that these ratios, below the 
QCD transition temperatures, are in reasonable agreement with a hadronic 
resonance gas. Above Tc the ratios rapidly drop and approach the high- 
temperature ideal-gas values. This suggests that, already for T > 1.5Te quark 
number and charge are predominantly carried by states with the quantum 
number of quarks. This leaves little room for the contribution of quark-quark 

^Note that this latter ratio could also include a contribution from quark-gluon bound 
states. Moreover, the observables discussed here arc insensitive to the presence of coloured 
quark-antiquark bound states, which could exist in the QGP [11]. 
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bound states. It is conceivable that these constraints would be even more 
stringent in 3-flavour QCD calculations, as the relevant states in the sQGP 
model are proportional to n'j while the relevant degrees of freedom in the 
(perturbative) high-temperature limit of QCD are only proportional to n/. 
A similar conclusion on the incompatibility of the sQGP model with lattice 
results on hadronic fluctuations has been drawn recently in [21]. 

It would, of course, be interesting to verify these results in experimental 
studies of event-by-event fluctuations [22, 23] of baryon number or charge. 
However, if indeed the medium created in a heavy-ion collision is strongly 
interacting and rapidly equilibrates there is little hope to find evidence for 
fluctuations originate from the QGP. At best, one may hope to study the 
quadratic and quartic fluctuations at freeze-out. Also this, however, would be 
interesting, as the observation of the importance of doubly charged hadrons 
leading to R4 2 — 1.5 could give further support to the validity of the reso- 
nance gas model. 
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